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ABSTRACT 

This report describes the research performed during the period March 
1986 - March 1987 under a joint research program between the Mechanical 
Engineering Department of the University of Maryland and the Center for Fire 
Research of the National Bureau of Standards. The research is conducted in 
the laboratories of the CFR by a Graduate Research Assistant of the ME 
Department under the joint supervision of Dr. Marino di Marzo (ME Dept.- 


UMCP) and of Dr. David D. Evans (CFR - NBS). 


The study of droplet evaporation on a high thermal conductivity semi- 
infinite solid is investigated as well as the solid thermal behavior both 
theoretically and experimentally. The behavior of a low thermal conductivity 
solid is the objective of experimental studies. A coupled analytical model 
for the solid and the evaporating liquid is being derived to predict the 
cooling effect induced by the evaporating droplet. The results of the 


analytical model will be available in the Summer 1987. 


In parallel with the single droplet study, a preliminary modelling 
effort is initiated to characterize the cooling effect of a multi-droplet 
array. Typical spray distributions are used to derive some information on 


overall cooling effect of a hot surface. 


A summary of the presentations and of the publications produced during 
the reporting period is also enclosed. A brief outline of the future research 
activities is given in the end of this report to provide a (prospective) of 


the reported results in the framework of the overall research program. 
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us INTRODUCTION 
eel Background 

The long term objective of the study of droplet-solid interactions is to 
obtain information applicable to the extinguishment of fire through a droplet 
array (e.g. spray). The solids of concern include porous materials, typical 


of fire applications. 


This report describes the research performed during the period March 
1986 - March 1987 under the joint research program conducted by the Center for 
Fire Research (NBS) and the Mechanical Engineering Department (University of 
Maryland) addressing the Transient Cooling of a Hot Surface by Droplet 


Evaporation. This joint research program was initiated in January 1985. 


The research described hereafter constitutes the first phase of an 
extensive research program aimed at developing accurate droplet cooling models 
of burning solid feet surfaces. Many studies have been performed to quantify 
the vaporization process of both single and multi droplet arrays impacting on 
hot surfaces. For the studies found in the published literature, the full 
span of the droplet vaporization processes are usually reported. These would 
include evaporation, nucleate boiling, film boiling and Leidenfrost 
transition. This present research is more limited in the span of vaporization 
process studied, being only concerned with the evaporation of single droplets 


on a hot surface. 


The results obtained in the first year of research are twofold. First, 


a simple mechanism of coupled heat conduction and evaporative mass transfer is 


identified as a plausible droplet air boundary condition. This information 
Ridews one to model the droplet evaporation in more complex situations where 
significant, two dimensional temperature gradients are present in both the 
droplet and the solid. Second, the formulation of a volume of influence in 
terms of a linear relationship with the evaporation time is the key step to 


derive a simple model for the cooling effect of a multi-droplet array. 


Lege Overview 
On the basis of these two important results, the research effort 


described in this report encompass several topics. 


The experimental results obtained in the first year of study show a 
strong dependence of the droplet initial configuration (on the solid surface) 
from the deposition technique. A serious attempt to standardize the droplet 
deposition procedure was undertaken in the first months of research. The 
details of the standardized droplet deposition procedure are given in the 


following sections. 


The results for solid of high thermal conductivity, obtained in the 
first grant period, compared well with similar studies available in the 
literature. However, the materials of concern in fire research exhibit low 
thermal conductivity. The information derived from the experiments on 
aluminum was applied to the study of the thermal behavior of glasslike 
materials. A first unsuccessful attempt to use regular glass was followed by 


a complete study of the thermal behavior of Macor. This material exhibits 


thermal properties similar to the regular glass but is resistant to the 


thermal shocks induced by the intense localized droplet cooling. 


The data collected in an extensive test matrix were then compiled in a 
spreadsheet format and plotted in a form suitable for the analysis. The 
ingredients necessary to develop a complex code for the prediction of the 
droplet and solid thermal behavior were identified. The governing equations 
were derived and a number of test cases were solved with the A.D.I. numerical 
"sethente? At that time also the Gauss-Seidel scheme was identified as an 


alternate numerical method to solve the governing equations. 


PSS Program Development 


In this section a brief synopsis of the various research activities is 


given in the time frame of the reporting period. 


Period March 1986 - May 1986 A new deposition device is designed, constructed 
and successfully tested. A new modular experimental set-up is also designed 
in order to provide simple interchange of different solid materials in the 
experimental rig. The A.D.I. numerical scheme is applied to a series of 


simple test cases. 


Period June 1986 - September 1986 The experimental rig is completed in its 


new configuration. Data for regular glass are collected but, subsequently to 
the mechanical failure of the glass block, the experiments on glass are 


terminated. An extensive data collection is performed on Macor. The A.D.1. 


numerical scheme is applied to the more complex case of cylindrical 


coordinates. 


Period October 1987 - December 1986 The data collected in the previous period 
are compiled in spreadsheet format and are plotted in a form suitable for the 
analysis. Note that during this period Mr. Meng successfully passed the 
written portion of his comprehensive examination at the University of 


Maryland. 


Period January 1987 - March 1987 The initial development of the code for the 
simulation of the thermal behavior of both the droplet Vand the solid is 
initiated. Simple boundary conditions are applied to test the accuracy of the 
numerical solution. In parallel with the work of Mr. Meng at the Center for 
Fire Research, Mr. Wang is hired to perform some preliminary computation on 


the modelling of the cooling effect of multi-droplet arrays. 


rah EXPERIMENTAL MEASUREMENTS 
2.1 Deposition Technique 

During the course of the aluminum data collection, it was observed that 
when a liquid droplet is deposited on a solid surface and spreads on it, its 
final configuration is affected by many factors and varies from one test to 
another. The volume of the droplet is measured by a micro dispenser which is 
used to gently deposit the droplet on the solid surface. It was found that 
the height of the end of the micro dispenser pipette from the solid surface 
hay the release height of the droplet) is a major factor affecting the 
final droplet shape. The angle at which the micro dispenser is held with 
respect to the normal to the surface and the cleanliness of the solid surface 


are other factors which influence the shape of the droplet. 


The evaporation time of a droplet that spreads more on the surface ae 
shorter than the evaporation time of a droplet that spreads less. Therefore, 
a deposition apparatus is necessary in order to strictly control the droplet 
release height and the position of the micro dispenser with respect to the 
surface. Furthermore, an effective technique to clean the surface must also 


be identified and an identical procedure must be followed for each test. 


Figure 1 shows the apparatus designed to deposit the droplet on the 
solid surface. The micro dispenser is mounted on a support which can be 
rotated clockwise of 90 degrees in the plane of the drawing. This enables the 
dispenser to be filled with water. Once the dispenser. is in vertical 


position, it can be locked in place. A motor driven cam actuates the plunger 
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of the micro dispenser. The cam is designed to execute a down stroke in half 
revolution followed by an up stroke which takes place in one quarter of a 
revolution. The last quarter of revolution leaves the plunger at its highest 
point. The motor speed can be controlled in order to obtain the most 


desirable deposition cycle. 


The height of the micro dispenser pipette with respect to the surface 
can be closely control by a fine pitch screw at the right end side of the 
apparatus arm. A substantial stand supports the arm and provides the 


necessary stability for the apparatus. 


The semi-infinite solid on which the droplet is deposited is constituted 
by a cylindrical block which has a diameter of 10.2 cm and height of 7.6 cm. 
The block is enclosed in a metal box shown in Fig. 2. The box is divided by 
an horizontal plate at the lower end of the cylindrical block. The upper 
portion of the box is filled with ceramic fiber insulation. The lower portion 
which contains the heater is empty and a number of large holes are drilled in 
its sides to insure proper ventilation of the heater. The bottom and the 
sides of the box are made of metal and both the top and intermediate partition 
are made of insulating, 1 cm thick, slabs. A hole of the same diameter of the 
cylindrical block is drilled in the top and intermediate plate. The 
cylindrical block is flush with the upper surface of the box and the heater- 


solid interface is at the intermediate plate level. 


The box contains also the thermocouples which are mounted on the solid 
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block. The box is mounted on an x-y positioned that allows the operator to 


exactly locate the droplet impingement point on the solid surface. 


ere Experimental Procedures 

In the previous year report [1], the properties of the water to be used 
in the experiments were discussed at length. The water is deionized to 
eliminate possible salt deposit on the evaporating surface. A complete 
degassing process is necessary because, otherwise, the droplet shape is 
altered by the release of the non condensible gasses during the evaporation 
transient. Furthermore, in the proximity of the nucleation temperature, the 
mon condensible gasses can quickly trigger nucleation, thus completely 


altering the nature of the droplet vaporization process. 


The surface is carefully cleaned first with 98 per cent ethyl alcohol 
and then with deionized water. Upon complete vaporization of the latter, the 
actual droplet is deposited. In order to deposit a droplet, the micro 
dispenser is filled with deionized and degassed water. The cam driven plunger 
pushes the water out of the micro dispenser pipette. Then the dispenser is 


rotated away from the surface. 


At the deposition time, a stopwatch is activated and a television camera 
monitors the droplet configuration. The evaporation time is recorded and 
information on the droplet shape parameter is deduced from the tv recording 


tape. 


se) 


The temperature gradient (in the vertical direction) across. the 
cylindrical block, prior to the droplet deposition, provides the necessary 
information to calculate the heat flux at the solid block upper surface. The 
convective and radiative components of the flux have not yet been determined 


due to the lack of information on the emissivity of the Macor surface. 


The droplet evaporation time, the shape parameter (e.g. the area of the 
wetted surface under the droplet), the time dependent history of the droplet 
shape and the overall heat transfer coefficient at the block upper surface are 


the data recorded during the various tests. 


fale Experiments With Glass 


The first low conductivity material selected for the experiments was 
common glass. The upper surface of the glass cylindrical block was polished. 
The bottom surface was covered with silicon heat sink compound to provide the 


necessary thermal contact with the heater. 


A few droplet were deposited on the block. A small crack due to the 
intense local thermal stresses was then observed. The crack propagated 


through the block slowly as more droplets were deposited on the surface. 
The clear failure of the common glass due to thermal stresses suggested 


more strict criteria for the selection of the low conductivity solid 


materials. Quartz and Macor were then selected. 
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2.4 Experiments With Macor 

The unsuccessful performance of common glass prompted the search for a 
material that could withstand the thermal shock induced by the localized 
droplet evaporative cooling. A machinable glass ceramic identified by the 
tradename: 'Macor’ seemed to have the required characteristics. Quartz was 


also found to satisfy the thermal stress requirements. 


The thermal and physical properties of Macor are illustrated in the Fig. 
3 through 6. The material was exposed to intense local thermal gradients and 
it performed according to the expectations. A cylindrical block with 
polished surfaces was used in the experiments. The dimensions of the block 
were 10.2 cm in diameter and 7.6 cm in height. The Macor block was mounted in 


the metal box as previously described for the glass. 


2.5 Data Collection 

Figure 7 illustrates the location of the various thermocouple on the 
Macor block. Two sets of three thermocouples are pointed at the top and 
bottom edges of the block equally spaced along the circumference. The 
temperature readings are used to infer the top and bottom surface 
temperatures. The heat flux across the block is then evaluated. In order to 
obtain information on the convective heat transfer coefficient it is necessary 
to evaluate the radiative component of the heat flux at the upper surface of 
the block. Measurements of the Macor emissivity are not available at this 
time. Figure 8 shows the behavior of the overall heat flux at the surface as 


a function of the surface temperature. Numerical values for the data points 


are also given in Table l. 
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Thermocouples 


* T;, To, T3 are temperature on the bottom 
surface 


ola tins ra) and Tg are temperature on the 
top surface 


¢ T7 is the ambient temperature 


e V; is proportional to the temperature difference 
between the top surface and ambient; and V> 
is proportional to the top and bottom surfaces 


FIGURE 7. Thermocouple Location on the Macor Block 
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TABLE ONE 
Overall Heat Transfer Coefficient 


Initial Solid Surface Overall Heat Transfer 


Temperature (°C) Coefficient (W/m* °C) 


The radius’ of the wetted area is measured from the enlarged tv picture 
of the droplet. An adequate scaling of the tv picture is performed for each 
test. The measured radius of the wetted area is ratioed to the radius of the 
spherical droplet of equivalent volume and this information, identified as the 
shape parameter §,, is presented in Fig. 9 for various droplet volume versus 
the initial surface temperature. The data base for the shape parameter 
includes 300 experimental data points. Each value shown in Fig. 9 and listed 
in Table 2 is the average of ten consistent experimental measurements. For 
the purpose of the analysis a constant value of the shape parameter of 1.42 is 
more than adequate since most of the experimental data are within a plus or 
minus five per cent band about this value. 
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TABLE TWO 


Shape Parameter B, 


 Aechenn. dot) Initial Solid Surface Temperature (°C) 
ee [ei [os [os [oo is 


1.34/1.44)1.42 


2371124571260 


The evaporation time for various droplet sizes and for various initial 
surface temperatures are listed in Table 3 and are presented in Fig. 10. The 
general trend of the evaporation time is in all similar to the trend of the 


aluminum data (see Fig. 11). Clearly, as expected, the cooling is more 


TABLE THREE 
Evaporation Time (seconds) 


Initial Solid Surface Temperature (°C) 


fname GoD] 96 foros [us Jus as [oo iss oo 
10 87 74 66 62 53 40 38 28 28 
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FIGURE 11. Droplet Evaporation Time Versus Initial Solid 


Surface Temperature (Aluminum) 
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intense and more localized. This means that the evaporation time for droplet 
deposited om Macor is longer than the time of evaporation of droplets 
deposited on aluminum given the same initial surface temperature. The reason 
for this result is evident if one considers that due to its lower thermal 
conductivity the Macor exhibits lower temperatures in the droplet proximity 
than the aluminum. In line with this consideration it is ico expected that 

nucleate boiling is suppressed for initial surface temperatures of the Macor 
in the order of 160 °C while, for droplets deposited on aluminum, nucleate 
boiling is observed at few degrees higher than the water saturation 
temperature. However, if one considers the contact temperature suggested by 
Seki [2] to be present under the droplet (at least at the initial contact 
time), it can be shown that, for an initial Macor surface temperature of 160 
°C and a liquid water initial temperature of about 25 °C, the contact 


temperature is very close to the water saturation temperature. 
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a THEORETICAL MODEL 
3.1 Theoretical Approach 

The experimental observation of the droplet evaporative phenomena can 
provide the necessary information to formulate a predictive model for the 
droplet evaporation. The validation of such a model is simple because both 
evaporation time and evaporation rate are easily deduced from the experimental 
records. More complex to validate and to correctly model is the behavior of 
the solid under the droplet. The energy balance equation governs the heat 
transfer in the solid; the boundary conditions to the energy equation are not 
clearly defined because a number of transport phenomena occurs at the solid- 
air interface (heat convection, radiation) and at the surface of the droplet 


(evaporative cooling, heat convection, radiation). 


A fundamental difference between high thermal conductivity solids and 
low thermal conductivity solids must be noted. The cooling is particularly 
intense at the solid-liquid boundary under the droplet. In a high 
conductivity solid a large portion of it around the droplet contributes to the 
heat supply for the droplet evaporation process because heat is transfered 
with relative ease under small temperature gradients. This is not so for low 
conductivity material where the limited heat transfer results in a more 
intense localized cooling in the immediate proximity of the droplet. Large 


temperature gradients are therefore expected in this case. 


The high conductivity of the solid allows one to assume that the 
temperature at the solid-liquid interface under the droplet is constant and 


equal to the contact temperature between two semi-infinite bodies as 
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demonstrated by Seki [2]. For the case of low conductivity materials this 
assumption is no longer valid and, in order to obtain a boundary condition 
under the droplet, the coupled solution of the energy equation for the liquid 


and the solid must be obtained. 


In summary, for the case of aluminum, a simple model to predict the 
droplet evaporation behavior is obtainable without solving the solid energy 
| equation whereas, for the Macor, the solid and the liquid energy equations 


must be solved simultaneously. 


A brief outline of the governing equation and their boundary conditions 
for the coupled case (Macor) is presented hereafter. Figure 12 represents the 
coordinate system used in the formulation of equations and boundary 


conditions. 


Solid energy equation 
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CNS Pine hn let) (3.4) 


) on) 


oT oT 
Ketiggic sBia5. (3.6) 
(T).= (D, (3.7) 


~k, SS =h, (T-1,) +0 € (TH Tr) | (3.8) 


@ y= f(r) (droplet-air interface) 


oT 
ab az h(t oe Oe vf 


ey oe 4 hy me 


Cy fiesta 


0.624 (2) ( 
A 


(3%9) 


where 0.624 is the air/water molecular weight ratio and f(r) is given by: 


2 
f(r) = R [ Wh. = i é 14 ie as A (3.10) 


and 


Gi 


abmrR (LEG CR PACU EEL bb wie ae 
3 Or 1/2. 173 
[ 4/pr— (1 + 16/89) 7717/7) (3.11) 
Recall that B is the shape parameter and is defined as the ratio of the 
radius of the wetted area (which is constant during the evaporation process) 
over the radius of the equivalent volume spherical droplet. Note that the 


volume of the droplet is changing as the evaporation takes place. 


Initial conditions 


(solid) 
oT 
_ Ko ao ha hy (T - T,) Sool eo 
(liquid) 
T= T, C1 5n 


T= ee (3.14) 
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3.2 Numerical Methods 
In order to solve the previous set of governing equations and boundary 


conditions two numerical schemes can be used. 


The first numerical scheme known as Alternate Directions Iteration 
(A.D.I.) can be applied separately to the liquid and to the solid. The 
relevant property of this scheme is the direct solution of the governing 
equation without iterations. This is due to the simple tridiagonal structure 
of the matrix to be inverted in order to obtain the solution of the algebraic 
system of equations obtained by discretizing the original governing equation 
(partial differential equation in time, radial and axial coordinates). The 
drawback of this scheme, however, is due to the iterative matching of the 
solid-liquid boundary conditions. Preliminary trials seem to indicate that in 
order to obtain a stable solution, one must resort to rather small time 
increments. Further study of the A.D.I. is necessary prior to obtain 


meaningful results with this numerical technique. 


The second numerical method is iterative but allows one to solve the 
liquid and the solid at once disregarding the internal boundaries (liquid- 
solid interface) and considering only the boundaries external to the liquid- 
solid system. This common numerical scheme is referred to as Gauss-Seidel. 
The advantage of the direct solution of the two regions is balanced by the 


less manageable iterative scheme. 


Both techniques are considered to solve the energy equation and a 


comparative evaluation of the results will be performed in the next grant 


period. 


Sis) 


3.4 


Nomenclature 

apex of the droplet (@ r = 0) 

specific heat at constant pressure 

air-vapor mass diffusivity 

equation of the liquid-air interface 
convective heat transfer coefficient 

thermal conductivity 

radius of the wetted area under the droplet 
radial coordinate tangent to the solid surface 


temperature 

initial solid surface temperature 
time 

vapor molar fraction 

air molar fraction 


coordinate normal to the solid surface in the liquid 


coordinate normal to the solid surface in the solid 


thermal diffusivity 

shape parameter 

emissivity of the solid surface 
latent heat of vaporization 
density 


Stefan-Boltzmann constant 
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Subscripts 


L 


Ss 


air free stream property 
liquid property 


solid property 


RP) 


4, MULTI-DROPLET MODEL 
4.1 Problem Statement 

The cooling effect induced by the evaporation of a single droplet on an 
aluminum surface has been measured and evaluated by diMarzo and Trehan [1]. 


Their study is the basis for the development of a multi-droplet model. 


The present study is aimed at developing a multi-droplet model on the 
basis of the single droplet information. Two questions must be addressed by 


the multi-droplet model. They are: 


1. How much water is sufficient for cooling the hot surface? 


2. What is the most efficient way of spraying water on the hot surface? 


Since the multi-droplet model incorporates experimental results of the 
single droplet experiment, the applicability of this model is confined within 
these experimental conditions. However, the concept and methodology employed 
in the model is not subjected to these restrictions. If a broader range of 
experimental data is available, the concept and methodology proposed in this 
study can be used to develop a similar multi-droplet model with wide 


applicability. 
In the following sections, a brief summary of the single droplet 


experiment and further analysis of the data are presented. Subsequently, a 


multi-droplet model is proposed and sample calculations are made. 


36 


ae Cooling of a Hot Surface by Single Droplets 


The cooling effect induced by the evaporation of a single droplet on a 
hot semi-infinite aluminum block was investigated by diMarzo and Trehan [1]. 
In this section, a brief summary of their results is presented. Also 
presented is a further analysis of their experimental data. This analysis 


leads directly to the development of a multi-droplet evaporation model. 


4.2.1 Summary of the Previous Results The study of diMarzo and Trehan [1] 
includes both the experimental and theoretical aspects. ra Shotts and 
the computer code simulation provide information on the cooling effect induced 
by an evaporating water droplet in a hot aluminum block simulating a semi- 
infinite body. The droplet is ‘gently’ deposited on the aluminum surface by a 
graduated hypodermic siring. The initial surface temperature ranges from 75 
°C to 100 °C; implying a full suppression of nucleate boiling. Thus, only 
evaporation takes place. The diameter of the spherical droplet prior to 
deposition on the surface varies from 0.8 mm to 4.5 mm. The detail of the 
experimental conditions can be found in Ref. [1]. In general, the following 


two main conclusions are the result of the study: 


ig The parameter §£,, proposed by Bonacina [4], is defined as the 
ratio of the droplet diameter after deposition (d) to that before 
deposition (D). Figure 13 illustrates its definition. Since the 
parameter B, characterizes the geometrical configuration ofthe droplet 
on the surface, it is an important parameter determining the droplet 
evaporating time and its cooling effect. It was found that if a droplet 


is gently deposited on an aluminum surface, the value of the parameter 


37 


Uo; IFUTZaq Jaqyauereg adeys “Et andi 


ews Aeids Burbuidwy e 
ne 6 ie jo|doip payisodap Ajjuay e 


anjen g jeoidAL 
q/p=g Jjejoweied 


38 


B, can be represented by a linear regression equation with an error bar 


of +/- 5% (Fig. 14). This equation can be written as 


B, = d/D = 0.6 + 0.009 * T. (4.1) 


where Le is the initial surface temperature (°C) 


ie The evaporation time was found to be an exponential function of 
the initial surface temperature and a power function of the droplet 
volume. A logarithmic regression analysis generates the following 
equation: 


(-0.07 T_) 
ie et ee ee ‘ (4.2) 


where +r is the droplet evaporation time (sec) and V is the droplet 
volume (ul). This equation yields values which are in the range of +/- 


20% of the experimental values. 


On the basis of the experimental data, a computer code was set up to 


calculate the thermal behavior of the semi-infinite solid in the neighborhood 


of an evaporating droplet. 


The concept of the volume of influence was proposed to quantify the 


cooling effect induced by an evaporating droplet. The volume of influence is 


defined as the volume of the region of the solid where the temperature drop 
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SHAPE PARAMETER By 


@ Experimental data 
— B,=0.009 T,+0.6 


70 75 80 85 90 95 100 105 
INITIAL SOLID SURFACE TEMPERATURE (°C) 


FIGURE 14. Shape Parameter Versus Initial 


Solid Surface Temperature 
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due to the droplet evaporation is at least 10% of the maximum possible 


temperature drop. This parameter is a function of several variables: 

1. The evaporation time 

2. The volume of the droplet 

3. The initial surface temperature 

4. The convective heat transfer coefficient 

It was found that a linear relationship exists between the maximum 
volume of influence (V/V) and the evaporation time (r). The relationship is 
represented by the following equation (See also Fig. 15): 


VV = 3+ 0.012 + (4.3) 


This linear relationship is fundamental for the developement of the multi- 


droplet model. 


4.2.2 Further Analysis On the basis of the volume of influence concept, two 
additional concepts can be defined in a similar way: the area and the radius 
of influence. The area and radius of influence are, respectively, the area 
and radius of the hot solid surface where the temperature drop due to a 


droplet evaporation is at least 10% of the maximum possible temperature drop. 


41 


NON-DIMENSIONAL 
MAXIMUM VOLUME OF INFLUENCE (Vj/V) 


50 100 150 200 250 
CALCULATED EVAPORATION TIME (s) 


FIGURE 15. Non-Dimensional Volume of Influence 
Versus Evaporation Time 
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An effort is made to correlate the maximum area and radius of influence with 


the droplet evaporation time. 


Figure 16 shows the maximum area of influence (A; ) as a function of the 
evaporation time for various droplet volumes and initial surface temperatures. 
It is found that no clear relationship exists between these two parameters. 

If the maximum area of influence is non-dimensionalized by the cross sectional 
area of droplet before deposition (Ap) > this non-dimensional parameter (A, /Ap) 
does not correlate with the evaporation time either (Fig. 17). However, if 
the maximum area of influence is non-dimensionalized by the area of droplet 
after depositing (A,)> the non-dimensional parameter (A, /Aq) correlates 
linearly with the evaporation time (Fig. 18). The linear regression of the 


data points can be expressed as: 
A,/Ag = 2.161 + 0.042 +r (4.4) 


Similarly, the non-dimensional maximum radius of influence (r,/T 4), also 


correlates linearly with the evaporation time (Fig. 19); that is: 
= 4 . 
r/t4 1.73 + 0.0075 +r (4.5) 


These two linear functions are important, since they are used directly 


to develop the multi-dropléet model. 
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4.3 Cooling of A Hot Surface by A Spray 

On the basis of the results of a single evaporating droplet, the cooling 
effect induced by a spray can be evaluated by taking into account the 
distribution of droplet sizes. In this study, a general model, which can be 
used for different types of droplet distribution functions and thus different 
types of sprays, is developed. By using this model, some sample calculations 
are also made to predict the minimum droplet mass per unit surface area 


necessary for the hot surface cooling. 


4.3.1 The Model In order to develop the general multi-droplet model, the 


following assumptions are needed: 
1. No coalescence of droplets takes place on the surface 


2. A suitable frequency function can describe the distribution of 
droplet sizes in the spray (the Rosin-Rammler [3,4], the Log-Probability 


[3] and the Nukiyama-Tanasawa [3,5] distributions). 


The droplet size distribution function can be written in several forms. 
In this study, the ‘volume distribution’ form is used; that is, dv/dD, where v 
is the volume fraction and D is the diameter of a droplet. Thus, the mass of 
droplets having diameter between D and D + AD that reach the surface between 


the time instant t and t + At is given by: 


dv 
A(Am) = (W At) (ap) AD (4.6) 
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where We is the mass flow rate of the droplet deposited on the hot solid 


surface. 


Therefore, the maximum area of influence produced by A(Am), Eq. (4.6) 


can be expressed by the following equation: 


A(Am) 
m 


A(AA, ) = F 


AA. (4.7) 


where m4 is the mass of a droplet and AA, is the maximum area of influence 
induced by a droplet. Equation (4.7) can be non-dimensionalized with respect 


to the total hot solid surface area A; that is: 


(4.8) 


For a single droplet, AA, /M 5 can be expressed in terms of the droplet 


diameter D, and the initial surface temperature Te It can be shown that: 


AA AA. A 
Say ae Nel 
7 Ars kh 
i em £2° 22 nr 3 
Pare (ira Da) / i Doe) (4.9) 


By substituting Eqs. (4.1) and (4.4) into the above equation, one can obtain 


the following result: 


AA, G2 LOL + 02042 er) CLS" (0.6 0.201 Te) 


5 (4.10) 


d Pa 
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Since r is given by Eq. (4.2), the above equation can be recasted as: 


AA, (-0.05 T_) 
See OL eee os DC a (4.11) 
are 


where c(T.) is a function Tos 12@a5 
iy pwn ey ES aha (60 + T Nahe (4212) 
s . s d i 


By combining Eqs. (4.6),(4.8) and (4.11), the following equation is 


obtained: 


AA. hee ipa Mie Ded 
i 0.092 s 
Se eC Ot one) a ora tae ie a0 (4.13) 


where We is the droplet mass flow rate in a spray. 


The total fraction of maximum area of influence produced by a spray of 
droplets having diameter ranging from Do to Da that reach the surface between 
the time t and t + At can be found by an integration of Eq. (4.13) over the 


droplet diameter range Ds to Ds: 


A. D (—0.05 T_) 
i m 0.092 cadv 
a ae AL) G(LS) Jp. 5 + Dive ] Gp oD 

fe dv dD 

= (Wy At) G(T.) [ 0.092 Jp a TET 

(<0 05"? 5.4) 
s m dv 
+e Jp Tape D dD } (4.14) 
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The two integration terms in Eq. (4.14) designate two different orders 

of mean diameters of a spray. In Ref. [3], the mean diameter Dae of order 
p + q is defined by: 

ces pede 

D x oh ay ae a aD] / ip pre ps (4.15) 

qP Do dD 

where An is the number fraction of droplets having diameter between D and D + 
AD. The mean diameters Dap can also be given in terms of Av; i.e., the 
volume fraction of droplets having diameters between D and D + AD. 


te Fis PPR tes 
f fee pq —- 3) < 3) “s 


Sak pe / i dDjee an 16) 


° 

Note that in many mass transfer and flow process [3], it is desirable to 
work only with mean diameters instead of the complete droplet distribution. 
Table 4 shows some of the mean diameters and some of their respective fields 
of application. These mean diameters can be readily calculated from the 


droplet size distribution equation. 


Consequently, according to the definition of the mean diameter, 
Eq. (4.12), the two integration terms in the Eq. (4.14) can be interpreted as 


and D respectively; that is, 


<a 
(D345) 43 


D D 
= fi 3 On m 2 dn 
27 Np Digi D1ie/ Up. D dD] 


dD 


D D 

m = m dv dadD_ Dn dv dv = 
e Sp p ddl / Jp ap sep sou / Wp ap dD ° ad 
(@] un 


ape 


TABLE FOUR 


Mean Diameters [3] (Cfr. Eq. 12) 


Linear Evaporation 


Surface Absorption 


Volume Hydrology 


Surface Dia. Adsorption 
Volume Dia. Diffusion 
Sauter Dia. Mass Transfer 


DeBrouckere Dia. {Combustion 


and, 


m dv m dv 
UR pee EUW A Taree cli 
Dn dv 
= lie Sy ale (4.18) 
Therefore, one can write 
AD Ga bieliley $e") 
1 OG092= = Ss 
ree (We At) c(T.) Ts + Di 3 € ] . (4.19) 
a2 
where os is commonly called the sauter mean diameter and a is referred 


to as the DeBrouckere diameter. 
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Note that these two mean diameters are related to the surface and volume 
of the droplets in a spray. The sauter mean diameter can be Paeeeoce cede 
the ratio of the total volume of the droplets in a spray to the droplets 
surface area; while the DeBrouckere diameter is the ratio of the square of 
the droplet surface area to the droplet volume. Consequently, if the two 
ratios, discussed above, for a single droplet are identical to those of a 
spray of droplets, the ratio A,/A calculated for a single droplet can be 


applied to the spray described by distribution function. 


Finally, the fraction of maximum area of influence is represented in 
terms of the droplet mass per unit area, (w At), the initial surface 
temperature, Ae and the mean diameters es and De The mean diameters can 


be calculated by using various droplet size distribution functions which 


represent different types of sprays. 


For the purpose of fire extinguishment, it is desirable to determine the 
minimum droplet mass necessary for cooling the hot surface, i.e., (wot) in’ 
In order to evaluate it, a criteria is needed to determine whether or not the 
hot surface is cooled. The proposed criteria is that the fraction of maximum 


area of influence is at least larger than 1.0, i.e., 


(4.20) 


IV 
re 


—  e2e where € 


Note that e« is a cooling index which indicates whether or not the hot 
surface is cooled. It is at least larger than 1.0, however, it should be 


determined by experiment. This cooling index is in principle similar to the 
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flooding index proposed by Rizza [6]. Furtherfore, the time interval At must 
be related to the average droplet cooling effect duration. The results of the 
solid thermal computation show that At can be set equal to twice the 


evaporation time of the average diameter droplet. 


By applying the above criteria to Eq. (4.19), (wo At) in can be 
calculated by the following equation: 


2. (—0.05 TS) ¢ 
(we At) vin Se { c(T.) (aD e a} (4.22) 


Hence, the minimum droplet mass per unit area for cooling the hot 
surface, (WoAt) in. is determined by the initial surface temperature, two 
different mean diameters and a cooling index. In the next section, some 
calculations are made to evaluate (wot) in for three droplet size 
distribution functions and for a range of initial surface temperatures, i.e., 


(Sp tOnLOU sso. 


4.3.2. Sample Calculation As noted in Section 4.3.1, (woAt) ints evaluated 
for the Rosin-Rammler, the Nukiyama-Tanasawa and the Log-Probability 
distribution functions. A brief description concerning these distribution 
functions is given. Subsequently, their application to the multi-droplet 


model is presented. 


An accurate understanding of droplet size distribution is fundamental 
for the analysis of the mass and heat transport. Because of its importance, 


numerous distribution functions have been proposed. The three functions 
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mentioned above represent three typical droplet size distributions of a spray 
which is considered as a system of liquid droplets in a continuous fluid 
phase. A critical examination of the droplet size distribution is referred to 


an excellent paper written by Mugelle and Evans [3]. 


A distribution function includes two essential concepts: a mean value 
and a measure of tne deviation from this mean. These two concepts are 
represented by the mean diameter D and the distribution parameter 6. By 
specifying proper D and 6, the distribution function can be used to predict 
the droplet size distribution in a spray. Subsequently, different order of 
mean diameters Dpcan be calculated from the distribution, for example, 


D35 and De which have been discussed in Section 4.3.1. 
The Rosin-Rammler (R-R) distribution function is generally given in the 


cumulative volume form: 


l-v 


wo 
~The (4.22) 
where 1 — v is the volume fraction of droplets having diameter larger than D. 
From Eq. (4.22), the volume distribution form can be obtained: 


So Tey) 


6 
cvreeor. A (4323) 


Paine rascé 


Lo] 


or in non-dimensional form: 
d 6-1 ak 
ss nee (4.24) 
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where 


z= D/D (4.25) 


The Nukiyama-Tanasawa distribution (N-T) function is an empirical 


equation. The volume distribution form is given as: 


dv ) pe/e 


) 
OTD UO ODO tag TO ge) 
dD [(6/6) Sooe 


D 


where b has the dimension of pan and [(6/6) is a gamma function of 6/6. If 
b is taken as D Cebu the N-T is represented by the following equation: 
(-6) 


oe, = re 
divin. tabbed 5 -(D/D) | 
aD 6750 (4.27) 


or in non-dimensional form: 


dv 0. ee pace (4.28) 


ADs Cer 6750 ee 


ol 


The Log-Probability (L-P) distribution function assumes a ‘normal 


distribution’ of y value; i.e., 


22 
EEL npn May ART BT (4.29) 


dy ae 


where y is defined as the natural logarithmic value of z; that is, 


y = ln z = 1n(D/D) (4.30) 
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In non-dimensional form is given as: 


5 ov ee 6 5 saan z)°] 


k Baan, d,dex) 
ee. 


Figures 20 through 22 illustrate the R-R, N-T and L-P distribution 
functions for three § values (i.e., 6 = 2, 3, 4). In these figures, the non- 
dimensional forms of the volume distributions are presented, Eqs. (24), 


(28) and (31). On the basis of these figures, two observations can be made: 


rigs When 6 is lower, the droplet size is more evenly distributed. On 
the other hand, when 6 is higher, the droplets in a spray tend to 
have the same size; that is, the width of the distribution curve 


is narrower. 


2: The most probable non-dimensional size (z), i.e., oe varies 


with 6 in each distribution. Table 5 lists Zn for the three 


distribution functions. 


The mean diameters of R-R, N-T & L-P can be evaluated by using Eq. 


(4.16). They are listed in Table 6 (see Ref. [2]). 
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TABLE FIVE 
The Most Probable z 


e 


TABLE SIX 


The Mean Diameter D 


(ptg-6) / (66 >) 
De 


Therefore, Dj, and D3 can be calculated by using the equations given 


in Table 6. They are presented in Table 7. 
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TABLE SEVEN 


abe Q(5/65) ] D Bil (T(6/6) ] 


ae 


In order to perform meaningful calculations by using the distribution 


functions, proper D and 6 have to be used. 


First, the criteria determining 6 is that the D/D range of the 
distribution function should be smaller than the single.droplet experimental 
range. In the experiment, the droplet diameter ranges from 0.8 to 4.5 mm; 


ive 


0.8<D<4.5 m (4.32) 
Ors De< 1) 7 (4.33) 
D = (0.8 + 4.5)/2 = 2.65 mm (4.34) 


Consequently, only the 6 = 4 curves for the three functions meet 
this criteria. Figure 23 illustrates three 6 = 4 curves. In order to 


compare them, the curves are shifted such that the most probable z, (D/D) 
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is equal to 1.0. The three curves show similar behavior. 


The determination of D is such that the most probable droplet diameter 
ae is the same for the three functions. The most probable diameter is 
taken as the average of the smallest (0.8 mm) and the largest (4.5 mm) droplet 


diameter in the experiment; i.e., 


Da (0.8 + 4.5)/2 = 2.65 mm (4.35) 


Therefore, D can be determined by the following equation: 
D my ala (4.36) 
where (D/D) ap £8 given in Table 5 for the three distribution functions. 


By using the values of D and the equations presented in Table 7, Dag 


and D3 are calculated. Table 8 gives D as well as De and be for the 


three distribution function with 6 = 4. 


TABLE EIGHT 
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Finally, by substituting the values of De Dao and De into Eq. 
(4.21), the minimum droplet mass per unit area needed for cooling the hot 
aluminum surface (wo 4t) 5, can be evaluated. Figure 24 illustrates three 
(wot) se ncurves for three distribution functions. Note that in this figure 


(wit) in is calculated by using the smallest cooling index; i.e., 


é = 1.0. On the basis of this figure, three conclusions can be drawn: 


pie When the initial surface temperature increases, (w At) vin 
increases. This implies more water is needed to cool the hotter 


surface. 


2h The three distribution functions give similar distribution. 


This is reasonable, since the three distribution curves (6 = 4) 
given in Figure 13 are similar though their functions are 


different. 
SF The predicted (w At). are in the range of 0.16 to 0.21 ee 


The conclusions discussed above are consistent and physically 
reasonable. However, the predicted (w At) sn needs to be verified by the 
multi-droplet experiment. It should be noted that the predicted value is 
obtained by assuming that « is equal to 1.0, i.e., € = 1.0. However, the 
cooling index « should be determined experimentally. The value of the time 
integral At should be evaluated as twice the evaporation time of the most 


probable diameter droplet (De oe 2.65 mm). 
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4.4 


A 


a 


Nomenclature 


total surface area 

area of droplet after deposit 

area of droplet before deposit 

maximum area of influence 

maximum area of influence induced by a droplet 
maximum area of influence induced by a droplet 
a parameter in the Nukiyama-Tanasawa Function 
a coefficient in Eq. (4.11) 

droplet diameter after deposit 


droplet diameter before deposit 
the largest droplet diameter in a spray 


the most probable droplet diameter in a spray 
the smallest droplet diameter in a spray 

a mean droplet diameter 

mean droplet diameter with an order of p+q 
Sauter mean diameter 

DeBrouckere diameter 

mass of droplets having diameter between D and 
D+AD that reach the surface between t and t+At 
Log-probability distribution function 

mass of a droplet 

number fraction of droplets having diameter 
between D and D+AD 

Nukiyama-Tanasawa distribution function 


an order index of mean diameter 
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q an order index of mean diameter 


ta maximum radius of a droplet after deposit 

tr maximum radius of influence 

R-R Rosin-Rammler distribution function 

t time 

At twice the evaporation time of the most probable diameter droplet 
T. initial surface temperature 

Av volume fraction of droplets having diameter 


between D and D+AD 


Vv initial volume of droplet 

V; maximum volume of influence 

We mass flow rate of droplets depositing on surface 
w mass flux of droplets depositing on surface 


(woAt) vin minimum droplet mass per unit surface 


for cooling hot surface 


y a variable 

fe non-dimensional parameter 
as the most probable 

€ a cooling index 

Pa density of the droplet 

re evaporation time of droplet 
‘ distribution parameter 

'h gamma function 
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High Thermal Conductivity Solid Surface", NBSIR 86-3384, 1986 


Mi di Marzo, D:. D. Evans, A. K. Trehan, "The Cooling Effect of a Single 


Droplet on a Hot Semi-Infinite Metal Body", NBSIR 87-3517, 1987 


evyrs Presentations 

Annual Conference of the Center for Fire Research (NBS Gaithersburg, October 
LvSG)": 

M. di Marzo, "The Cooling Effect of a Single Evaporating Droplet on a Hot 


Semi-Infinite Metal Body" 


Eastern Session of the Combustion Institute (Puerto Rico, December 1986): 


M. di Marzo, D. D. Evans, "The Cooling Effect of a Single Evaporating Droplet 
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High Conductivity Solid Surface", Accepted for presentation at the 24th 


ASME/AIChE/ANS National Heat Transfer Conference, Pittsburgh, 1987 
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6. SUMMARY AND PERSPECTIVES 
6.1 Low Conductivity Solid 

The experimental results shown in the previous section are in all 
similar to the results obtained for aluminum. A notable difference is in the 
behavior of the shape parameter B,: The data for aluminum indicated an upward 
trend of this parameter with rising temperatures, while, for the Macor, the 


data seem to show a constant value of about 1.42 (+/— 0.09). 


The results of the code for the coupled liquid solid boundary will give 
some important indication concerning the radial temperature gradients at the 
liquid-solid interface. Furthermore, the code will also validate the 
assumption of uniform temperature under the droplet used for the high thermal 


conductivity solid (see also Seki [2]). 


6.2 Radiation 

Upon successfully demonstrate the capability of the code to predict the 
droplet-solid interaction, the more complex situation where the major heat 
source is radiant in nature can be addressed. In a typical fire application, 
the heat is radiated from the flames to the exposed solid surfaces. The 
droplet impinging on such surfaces will receive radiant heat from above as 
well as heat from below conducted through the solid. In order to analyze this 
situation, the solid-liquid model will be upgraded to incorporate the 
radiative heat transfer component. The numerical solution to this complex 
problem will be validated with adequate experimental measures. Both the 
radiative model and experiments will be performed in the second portion of the 


next grant period (1987-88). The results of this investigation will be 


71 


hopefully collapsed in simple area (volume) of influence relationships similar 


to those derived for the simple aluminum case. 


heer! Spray Cooling 


The long term objective of this research program is to obtain a reliable 
model for solid fuel fire extinguishment. The results of the radiative model 
and experimentation will provide the basic information to derive a prediction 
of the cooling effect of a multi-droplet array. The methodology to extend the 
single droplet results to a droplet array is outlined. Validation of these 
predictions could be then derived from spray cooling measurements in actual 


fire environment. 
The correct prediction of the cooling induced by a spray on solid 


surfaces exposed to fire will constitute a preliminary milestone in the study 


of the complex problem of solid fuel fire suppression. 


72 


[3] 


[4] 


[5] 
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